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Abstract 

We treat the problem of a Michelson interferometer in the field of 
a plane gravitational wave in the framework of general relativity. The 
arms of the interferometer are regarded as the world lines of the light 
beams, whose motion is determined by the Hamilton-Jacobi equation 
for a massless particle. In the case of a weak monochromatic wave we 
find that the formula for the delay of a light beam agrees with the result 
obtained by solving the linearized coupled Einstein-Maxwell equations. 

We also calculate this delay in the next (quadratic) approximation. 

1 Introduction 

The subject of this work is to investigate the behavior of a Michelson inter¬ 
ferometer in the presence of a plane gravitational wave. Sources of gravita¬ 
tional waves include the collapse of massive stars in supernova explosions, 
and the coalescence of stars in binary systems containing a neutron star 
and/or a black hole. The large distances from these objects allow us to 
assume that the gravitational waves detected on the Earth may be regarded 
as plane waves. 

Gravitational waves were predicted by Einstein in his general theory 
of relativity PEI as perturbations of the spacetime metric that propagate 
at the speed of light. Plane gravitational waves were considered first by 
Rosen jHj , Taub P , and McVittie P , and their work seemed to indicate that 
such waves may not physically exist. However, Bondi, Pirani, and Robinson 
PPP proved that this conclusion was wrong, and that these waves not 
only exist but also carry energy. Exact solutions of the Einstein equations. 
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which describe plane gravitational waves, have also been developed in 
nni. In 1975, Hulse and Taylor found indirect proof for the existence of 
gravitational waves by observing the pulsar PSR 1913+16 HU. They noticed 
that the orbital parameters of this pulsar were changing according to a 
model assuming that such a system radiates gravitational waves urns]. 
This discovery was awarded with the Nobel Prize in physics in 1993. 

Research on the possibility of detecting gravitational waves was started 
in the 1960s by Weber HU, who constructed the first resonant-mass antenna. 
His results claiming the discovery of gravitational waves m could not be 
confirmed by other groups who built similar detectors. In the 1970s, another 
method of detection was suggested using laser interferometry HEllIZlIIHlIin]. 
Current and planned ground-based laser interferometric detectors include 
AIGO (Australia) [201, GEO 600 (Germany/UK) [2U, LIGO (USA) [22], 
TAMA 300 (Japan) [2Sj, and VIRGO (Prance/Italy) HU- The first search 
for signals from neutron stars in LIGO and GEO 600 data did not find direct 
evidence for gravitational waves but gave upper limits on the strength 
of periodic gravitational radiation HU- In 2011, ESA and NASA plan to 
launch LISA a space project that will be able to make observations in 
a low-frequency band that is not achievable by ground-based observatories. 
A bibliography on gravitational wave theory and experiment through 1999 
is listed in EHj. 

Since the search for gravitational waves is mainly based on interferome¬ 
try, we consider a Michelson interferometer in the presence of a plane gravi¬ 
tational wave. The effect of such a wave on a light beam, using the linearized 
coupled Einstein-Maxwell equations, was studied in EiEni. However, these 
equations in the quadratic approximation may lead to quite lengthy expres¬ 
sions. Instead, we can treat the light beam as photons whose motion is 
determined by the geodesic equations for a massless particle. 

The equations of motion for a particle in the metric of a strong grav¬ 
itational wave m were derived by Ebner m, and all seven integrals of 
motion can be found in the literature, e.g. m A simple model of a Michel¬ 
son interferometer in the field of such a wave, based on the above equa¬ 
tions of motion, was proposed by Bazahski m- Here, we generalize this 
model to the second (quadratic) approximation. For simplicity, we consider 
a monochromatic plane gravitational wave whose direction is parallel to one 
arm of the interferometer and perpendicular to the other. Since the energy 
of gravitational waves contains the factor ^, the amplitude for such waves is 
very small. Therefore, second-order effects in gravitational-wave detection 
are irrelevant experimentally, although it might be interesting to explore 
them. 
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This paper is organized as follows. In section 2 we introduce a plane 
gravitational wave and review the equations of motion for a test particle. In 
section 3 we study the behavior of a Michelson interferometer under the ac¬ 
tion of this wave in the first (linear) approximation. In section 4 we proceed 
to the second approximation. The results are summarized in section 5. For 
the metric tensor and the curvature tensor we use the conventions of CSl. 

2 Motion of a particle in the field of a plane grav¬ 
itational wave 

In this section we derive the motion of a test particle in the presence of 
a plane gravitational wave isn ES] using the relativistic Hamilton-Jacobi 
equation. Let us consider the metric tensor whose components are functions 
of only one variable u, gab = gah{u)- In this case, it can be shown that 
the spacetime interval in vacuum is of the form 

ds^ = dudv + gab{u)dx°‘dx^, ( 1 ) 

where the letters a, b refer to the coordinates 2,3. For a plane wave propa¬ 
gating in the direction of the x axis we have u = ct — x and v = ct + x. If 
the interval is given by equation m then all components of the Ricci tensor 
identically vanish except Ruu = where = Qab, and the dot 

denotes differentiation with respect to u. Thus, the Einstein field equations 
in the case of a plane gravitational wave reduce to one equation, Ruu = 0. 

The Hamilton-Jacobi equation for a particle with mass m moving in the 
gravitational field described by the metric in 0 is given by 

, dS dS dS dS 

4-h Q ° -= 

du dv dx^ 

where S = S{u, v, x“) is a principal function. We seek a solution of ((2j in the 
form S = S{u) + \Bv + where B = pQ + = const and poa = const. 

For a timelike geodesic line we have B > 0. We obtain the principal function 
given by 

S = ^[{rn^c^u - G“^(M,0)poaPob)] + ^Bv + poax'", (3) 

where we define 

G^\p,q)= r%'^\u)du. (4) 

Jq 
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The components of four-momentum are given by the derivatives of a 
principal function, pi = The world line of a particle is determined by 
Jacobi’s theorem, = /3^, where are the constants in the principal 
function 0, and are new constants. Introducing a finite interval s as 
a parameter such that yields u = — so) + ^o- The 

equations of motion for a massive particle in the field of a plane gravitational 
wave in terms of s are ^ 


ct{s) 


ct{so) + - 


f B mc\ 
\mc ^ B ) 


{s - So) 


1 

x{s) = 


cab 


B 


me 


{s - so),ct(so) 


x(so) + 


1 




- x{so)jpOaPOb, 
{s - So) 


~ So),Ct(so) - xisoi^POaPOb, 

x“(s) = x“(so) + —(s - so),ct{so) - x(so) jpofe- 

B \mc J 


( 5 ) 


In a similar manner, we can determine the equations of null geodesic 
lines in the presence of a plane gravitational wave. The Hamilton-Jacobi 
equation in this case is equivalent to the eikonal equation, 

and its solution is of the form 'I' = T(u) J- \Cv + koaX^", where C = ko + 
kx = const and koa = const. Consequently we find 

T = 0)koakob + ^Cv + feoaX^, (7) 

with defined in Q. The components of the four-dimensional wave 
vector are given by the derivatives of the eikonal, ki = A world line for 
a massless particle is again determined from Jacobi’s theorem. 

Let us define the affine parameter vr such that = 5 “^^^ from which 
we find u = C{'k — ttq) + uq. As a result, we obtain the null geodesic lines: 

ct(7r) = ct(7ro) -h ^C'(7r - vro) 

^In agreement with m 
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- TTo), ct(7ro) - x(7ro))koakob, 
x(7r) = x(7ro) - ^C(7r - tto) 

- vro), ct(7ro) - x(7ro))koakob, 

x“(7r) = x“(7ro) + ^G“'’(C'(7r - tto), ci(7ro) - x{To))kob, (8) 

where G°‘’^{p,q) were defined in (jl]). Formulae (jSl) are satisfied for all null 
geodesics, except the curve for which C = 0. This condition represents a 
light beam moving parallel to a gravitational wave. In this case, ka = 0 for 
an arbitrary vr, and the solution is 

ct = A{'K — TTq) + CtQ, X = A{'K — TTo) + Xo, x“ = Xq, (9) 

where A is a constant. 

3 A Michelson interferometer in the presence of a 
plane gravitational wave 

Using the equations of motion for a light beam in the field of a plane grav¬ 
itational wave, Q and ®, we will show how to build a simple model of a 
Michelson interferometer !MI. Let us consider three observers at rest: 

1. ct = s — So, X = 0 , y = 0 , z = 0 , 

2 . ct = s — So, x = ljj, y = 0 , z = 0 , 

3. ct = s — So, X = 0, y = l±, z = 0. (10) 

Let a light beam leave the world point (0,0,0,0) in the direction of a plane 

wave, and arrive at the point (si — so,/||,0,0). In this case equations Q 
lead to 

Si - So = ^(ffi - ffo), /|| = A{^i - 7fo), (11) 

which gives si — sq = l\\- Then, let this light beam leave the point (si — 
So, ^ 11 , 0 , 0 ) and arrive at the point (s 2 — so, 0 , 0 , 0 ) (return to the initial space 
point). Consequently, the equations of motion Q are 

S2 - Si = ^C(7r2 - TTi) - ^G“'’(C'(s 2 - Si), Si - /||)/sOa^Ob, 

“^11 = - 2 ^(s 2 - Si) - -^^G‘^^{C{s2 - Si), Si - l\\)koakob, 

0 = ^G<^\C{s 2 - Si), Si - l\i)kob. (12) 
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Therefore S 2 — sq = 2Z||. 

Next, let us consider a light beam leaving the point (0,0,0,0) in the 
direction perpendicular to that of the gravitational wave, and arriving at 
(si — so,0,/_L,0). In this case equations (|H1) give 


Si — •So = 2^^^^ ~ ~ -^^G'i^koakob, l± = -^G'^kob, 

0 = -\g{tt - TTo) - 0 = ^G^Jikob, (13) 


where 
for si: 


G'^’^{C{7ri — ttq),so) = g°'^{u)du. We obtain the equation 


Si — So = — 


llG^J^ 


(14) 


where - {G^J^f. 

Then, let the same light beam go from the point (si,0,/_l, 0) to the 
initial space point (s 2 ) 0,0,0). The equations of motion are (with different 
constants of motion) 


S2 - Si = ^<^(712 - TTi) - ^^Gt.kQakob, -l± = 

0 = -^^(712 - TTi) - A:oaA:ob, 0 = ^G^kob, (15) 


where G^ 
obtain 


G“^((7(712 — 7ri),si) = g'^^{u)du. In a similar manner we 


S2 - Si = - 




(16) 


where = G^^G^J^ - (G^)^. 

A weak plane gravitational wave in metric o can be described by the 
two-dimensional metric tensor 


-l + f{u) h{u) 
h{u) -1 - f(u) 


(17) 


where f(u) and h(u} are small quantities The two possible polarizations 
correspond to setting either quantity to zero. Let us consider a polarization 
such that h(u) = 0 and denote it by /i+. From this we find A^ = (si — 

^If we neglect terms quadratic and smaller in / and h, then the metric tensor m 
identically satisfies the field equation Ruu = 0. 
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■so)^ + 0{A'^), where A is the amplitude of the gravitational wave. We can 
linearize equation m, omitting terms quadratic and smaller in A: 

^0 = so) + f{u)du^. (18) 

The analogous expression can be derived for equation 

Let us assume that a plane gravitational wave is monochromatic with 
wavelength A, f{u) = Acos^^. We seek the solution of equation (1181) in 
the form si — sq = where As_> is a small quantity on the order 

of A. Keeping only linear terms in As^ we find 


AX . 7rl± 27r(so + 

As^ = - sm——cos -— 

27r A A 

Similarly, for the returning beam we have S 2 — si = /^ + As^, where 



(19) 


As 


AX Trl± 27r(so + 3 /_l/2) 

- sin——cos -;-. 

27r A A 


( 20 ) 


where we used si — sq = li_ since this term is on the order of A. Finally, 
the total delay of the light beam moving perpendicular to the gravitational 
wave, relative to the parallel beam, is given by 


As = As^ + As 


AX 27r(so + ^_L) 

- sin—;—cos -;-. 

27r A A 


( 21 ) 


The other polarization (denoted by /ix) is described by the metric tensor 


—1 h{u) 
h{u) —1 


( 22 ) 


In this case, equation m reduces to si — sq = whereas for the returning 
beam we find S 2 — For this polarization there is no delay. 

Equation m is what experimentalists use to detect gravitational waves 
[2H1IS1- Its relation to incoming gravitational waves with arbitrary prop¬ 
agation directions and polarization states can be found in the literature, 
e.g. 122]. This equation agrees with the corresponding expression in isni if 
we set the departure moment sq = 0 and use the relation As = ^A(p, where 
LV is the light beam frequency and A4> denotes the phase shift. Therefore, 
treating light as massless particles in the field of a weak gravitational wave 
and solving the linearized coupled Einstein-Maxwell equations lead to the 
same result. 
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The delay of a light beam depends on the quantity sq, which corresponds 
to the moment when the beam leaves the start point (the time dependence 
of the gravitational field is fixed by the form of f{u)). If the wavelength 
of a monochromatic gravitational wave is small compared to the size of the 
interferometer, we can average this delay over sq: (As) = j As(so)dso- 
The physical reason for averaging arises from the fact that we do not know 
when the front of a gravitational wave hits the point of the departure of the 
beam. However, as for each periodic function, such an expression vanishes 
and is of no interest. 


4 Light beam delay in the second approximation 


In this section we will generalize the results of the preceding section and 
consider quadratic terms. Let the two-dimensional metric tensor gab of a 
weak gravitational wave be given by 


9ab — 


hw -1 + ^(1)]’ 


(23) 


where are small quantities on the same order. In the linear 

approximation where A is the amplitude of such a wave 

on the order of Therefore, without loss of generality we may assume 
= 5 ^^^- The field equation is obtained from Run = 0: 

1(/(1)2 + 2/(i)/(i) + + /,(1)2 + = 0, (24) 

which gives < 7 ^^^ as a function of and 

If we consider a polarization /i+ with = 0 and put si — sq = L in 
terms on the order of A^, then equation (1141) becomes 

isi-sof-l‘i{si-so)+lif + f f^^'^du'] =0, (25) 

J So Sq / 

Note that this equation does not contain the quantity . Let us apply the 
above results to the case of a monochromatic plane wave, = Acos^^. 
For such a wave, equation leads to — 1), and this 

quantity clearly is not a periodic function of n We seek the solution in 

®We could assume but then equation. 12511 would contain the nonpe¬ 

riodic function Thus, it would be impossible to average the light beam delay over sq 
(the moment when the beam leaves the start point) independently of the initial conditions. 
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the form si — sq = Z_l + As^, where As^ is again a small quantity on the 
same order as A. Thus, up to quadratic terms in As^ we obtain 


. . 7rl± TT . A^\ . 27rl± 2 tt . , , 

As^ =- sm——cos—{2so + - sin —-—cos—(2so + t±) 

27r A A Ihvr A A 


+■ 


AH^ 

4 

A^A 


+ 


A^A^ 

Stt^I I 


7rlj_ TT 

sin——cos—(2so + i±) 
A A 


7r/j_ TT , , 27r . , 

szn——cos—(2so + l±)cos—(so + l±). 
47 r A A A 


(26) 


The above expression contains the quadratic corrections to equation (HU). 
The mean value of this quantity over the moment of the departure of the 
beam sQ) = j As^(so)dso, is now different from zero: 


(As^) 


^2a2 

IGvr^T 


sin^ 



A^l^ 

4 


27r/j 

+ -r^sin- 


IGvr 


(27) 


We remind that the moment sq is relative to the instant when a front of 
the gravitational wave in question hits the departure point, and we average 
with respect to sq since we do not know the value of this quantity. 

For the returning beam we must replace si with S 2 and sq with si (since 
the beam departs from the reflection point right after it arrives there), which 
gives S 2 —Si = Z_L +As^, where As^ is described by equation (OHl) in which 
So is replaced with sq + l± + As^. The result is 


IA \ /A \ 27 r/_L 

(As^) = (As^)- :—sin ——sin —;—. 

47r A A 

Therefore, the total average delay of the light beam is given by 

'2'2 A'^L 


,A N ^ ^ • 2 ' 

— ■ 

A^A . oTtI± . 2ttI± 

- sin ——sin —;—. 

dvr A A 


212 A 2Trl 
+ —sin 


Stt 


A 


(28) 


(29) 


We observe that (As) is negative for all values of A. If Zy = l±, then the 
light beam moving perpendicular to the direction of a gravitational wave 
lags (on average with respect to the departure moment sq) behind the one 
which moves parallel, provided both beams left the light source at the same 
time. The dependence of the quantity D = on the variable / = is 
shown in figure ^ 
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Figure 1: The normalized average light beam delay D as a function of the 
normalized wave frequency / for the polarization 


Now we consider the other polarization hx for which = 0. Equa¬ 
tion (1251) is replaced by 

(si — so)^ —/l(si — So) + — l±( f h^^^du\ =0. (30) 

J so \Jso / 

For a monochromatic plane wave, = Acos^^, instead of equation (HH) 


we find 


As^ = 


f . 'kI±_ tt 


. 27r/x 


sin—^cos—{2so + l±)j - ^ 


A^h 


27r 

-sm^—cos—{2so + l±). 
Ibvr A A 


(31) 


The returning beam again satisfies formula (P|) . Averaging over the depar¬ 
ture moment sq gives 

A^A^ . . 2^^-L • 27r;_L ^ ^ 

' ^ 2tt^I I a 2 47r a a ^ ^ 


We see that (As) is negative for all values of the length of a gravitational 
wave, and tends asymptotically to zero as A —> oo. In figure we show the 
function D{f). Note that for both polarizations D —s- 1/2 as A ^ 0. 

From figures Hand [21 we conclude that a one-to-one relation between the 
average light beam delay and the wavelength exists only for waves longer 
than ~ 7r/_i_. Moreover, the precision of the measurement of the delay is 
proportional to the length of the arms of the interferometer. Therefore, if 
we try to increase the precision by increasing this length, we end up with 
decrease of the range of waves measurable in this way. 
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Figure 2; The normalized average light beam delay D as a function of the 
normalized wave frequency / for the polarization hx ■ 

5 Summary 

In this work we showed that the time delay of a light beam in a Michelson 
interferometer, obtained by treating light as massless particles moving on 
geodesics, agrees with the corresponding solution of the linearized coupled 
Einstein-Maxwell equations. We also calculated this delay in the second 
(quadratic) approximation. We considered the simple case of a plane wave 
with the direction parallel to one arm of the interferometer, and treated the 
two polarizations separately. However in reality, one should expect a broad 
spectrum of frequencies and a combination of both polarizations. 

Averaging the light beam delay As over the moment of the departure of 
the beam sq would make sense if either wavelength A is much smaller than 
the arm length l± or we measure (As) many times. Since we excluded the 
former, this method would not work in the case of supernova explosions that 
radiate significantly for a few seconds only. Instead, the presented model 
could be used (in principle) for continuous sources of gravitational waves 
such as stellar binaries containing a neutron star. As we stated in section 1, 
second-order effects in gravitational-wave detection are extremely small and 
thus irrelevant experimentally. However, they might be of some academic 
interest. 
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